Numerical Normal Forms for Codim 2 Bifurcations of Fixed Points with at Most Two Critical Eigenvalues

In this paper we derive explicit formulas for the normal form coefficients to verify the nondegeneracy of eight codimension two bifurcations of fixed points with one or two critical eigenvalues. These include all strong resonances, as well as the degenerate flip and Neimark-Sacker bifurcations. Applying our results to n-dimensional maps, one avoids the computation of the center manifold, but one can directly evaluate the critical normal form coefficients in the original basis. The formulas remain valid also for n = 2 and allow one to avoid the transformation of the linear part of the map into Jordan form. The developed techniques are tested on two examples: (1) a three-dimensional map appearing in adaptive control; (2) a periodically forced epidemic model. We compute numerically the critical normal form coefficients for several codim 2 bifurcations occuring in these models. To compute the required derivatives of the Poincaré map for the epidemic model, the automatic differentiation package adol-c is used.